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Abstract. After recent improvements to the Pavia model of two-nucleon knockout from 16 O with elec- 
tromagnetic probes the calculated cross sections are compared to experimental data from such reactions. 
Comparison with data from a measurement of the ie O(e,e'pn) reaction show much better agreement be- 
tween experiment and theory than was previously observed. In a comparison with recent data from a 
measurement of the le O(7,pn) reaction the model over-predicts the measured cross section at low missing 
momentum. 

PACS. 21.30.Fe Forces in hadronic systems and effective interactions - 21.60.-n Nuclear structure models 
and methods - 25.20.Lj Photoproduction reactions - 25.30.Fj Inelastic electron scattering to continuum 



1 Introduction 

A major quest of nuclear physics is to understand how 
the properties of nuclei arise from the underlying nucleon- 
nucleon (NN) interaction. A useful starting point is given 
by the independent particle models (IPM), in which pro- 
tons and neutrons move freely in a common mean field. If 
■ one accounts for spin-orbit effects and the average effects 
of the tensor interaction [T], this approach explains the 
shell ordering of most stable and dripline isotopes. How- 
ever, this picture cannot describe other basic observations, 
such as the strong fragmentation of nuclear spectra and 
the corresponding quenching observed for absolute spec- 
troscopic factors 2 . The failure of the IPM arises from 
the correlated behaviour between nucleons, which, at short 
inter- nucleon separations, is characterised by a strong re- 
pulsion and, at intermediate to long range separations, by 
an attractive interaction dominated by complicated tensor 
and spin-orbit terms. Thus, to understand nuclear struc- 
ture a careful study of this correlated behaviour is vital 

HI- 

A direct method to study NN-correlations is by the 
use of two-nucleon knockout reactions with an electromag- 
netic probe [5] . Proton- proton and proton- neutron knock- 
out reactions can act to probe the short range and tensor 
components of the NN-interaction, respectively. Real and 
virtual photons provide different and complementary in- 



formation on the reaction process. Real photons are only 
sensitive to transverse components of the interaction while 
virtual photons are sensitive to both the transverse and 
longitudinal components. 

Electromagnetically induced two-nucleon knockout re- 
actions are driven by several processes. The coupling of 
the (real or virtual) photon to either nucleon of a corre- 
lated pair via one-body hadronic currents can lead to the 
ejection of both nucleons from the nucleus. Interaction of 
the photon with two-body hadronic currents such as me- 
son exchange currents (MEC) or isobar currents (IC) also 
contributes to the cross section. In addition final state in- 
teractions (FSI) between the two ejected nucleons and the 
recoil nucleus need to be taken into account. The relative 
importance of these different processes depends on the re- 
action type and kinematics. 

The 16 O nucleus is of particular interest regarding the 
study of correlated behaviour. Various theoretical mod- 
els exist which attempt to describe these reactions [HE 
H] and there have been numerous measurements of two- 
nucleon knockout reactions using both real and virtual 
photons pfTUl[TTlll21ll31ll41ll5j . In the results from the 
first measurement of the 16 0(e,e'pn) 14 N reaction were re- 
ported and compared with theory. Theoretical predictions 
for this reaction have been obtained in [15] by combining 
the self-consistent Green's functions method for correla- 
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tions and the Pavia model for the reaction mechanism. 
These calculations suggest a strong sensitivity of the cross 
sections to tensor correlations. However, they were unable 
to reproduce the shape or the magnitude of the data [9]. 
These discrepancies sparked further developments to im- 
prove the reaction model with respect to the treatment 
of FSI [P7lfT8]. of the two-body currents [19], and of the 
centre-of-mass (CM) effects in the electromagnetic cur- 
rent operator [6ll20] . This paper presents a new compar- 
ison between experimental data and recent calculations 
[BJ, and shows that CM effects resolve the discrepancy 
found in [9]. We also show calculated cross sections for 
the 16 0(7,pn) 14 N reaction and compare these to a recent 
measurement of this reaction [21] . 

2 Theoretical calculations 

The cross section of a reaction induced by a real or virtual 
photon, with momentum q, where two nucleons are ejected 
from a nucleus can be written in terms of the transition 
matrix elements of the nuclear current operator between 
initial and final nuclear states. Bilinear products of these 
matrix elements give the components of the hadron ten- 
sor and therefore the cross section [5]. For an exclusive 
process, where the residual nucleus is left in a discrete 
eigenstate of its Hamiltonian, and under the assumption 
of a direct knock-out mechanism, the transition matrix 
elements contain three main ingredients: the two-nucleon 
overlap function between the ground state of the target 
and the final state of the residual nucleus, the nuclear cur- 
rent, and the two-nucleon scattering wave function [22] . 

The two-nucleon overlap function (TOF) contains in- 
formation on nuclear structure and correlations. In [B] dif- 
ferent treatments of correlations are compared, and pro- 
duce dramatic differences both in the shape and in the 
magnitude of the proton-neutron emission cross sections. 
In particular, a crucial role is played by tensor correla- 
tions. In the most refined approach, the TOF is obtained 
from a self-consistent calculation of the two-hole Green's 
function. In this case, the coupling of nucleons and collec- 
tive excitations of the system is calculated microscopically 
from realistic NN forces. This is done employing the Fad- 
deev random phase approximation (FRPA) method dis- 
cussed in [2Ull2"l] . The long-range part of tensor correla- 
tions is also included explicitly. The TOF has been calcu- 
lated in [16] by partitioning the Hilbert space. Long-range 
correlations are evaluated using FRPA and the Bonn-C 
NN-potential [25[[26] in an appropriate harmonic oscilla- 
tor basis. The effects of short-range correlations, due to 
the central and tensor part at high momenta, lie outside 
this space. Thus they were added by computing the ap- 
propriate defect functions. 

The nuclear current is the sum of a one-body and a 
two-body contribution. The one-body current includes the 
longitudinal charge term and the transverse convective 
and spin currents. The two-body current is derived from 
a non relativistic reduction of the lowest-order Feynman 
diagrams with one-pion exchange and includes terms cor- 
responding to the 7r-seagull and pion-in- flight diagrams, 



and to the diagrams with intermediate Z\-isobar config- 
urations. Details of the nuclear current components can 
be found in [19,27,28J. In comparison with the previous 
calculations of [9] , the treatment of the two-body current 
has been improved using a more realistic regularised ap- 
proach, which is consistent with the description of elastic 
NN-scattering data [BJ 

The two-nucleon scattering wave function contains the 
interaction of each one of the two outgoing nucleons with 
the residual nucleus, described in the model by an opti- 
cal potential, as well as the mutual interaction of the two 
ejected nucleons (NN-FSI). In a simpler approach (DW) 
only the contribution of FSI due to the optical poten- 
tial is included, and the scattering state is written as the 
product of two uncoupled single particle distorted wave 
functions, eigenfunctions of a complex phenomenological 
optical potential which contains a central, a Coulomb, and 
a spin-orbit term. In the more complete approach (DW- 
NN) the contribution of NN-FSI is also included within 
the perturbative approach reported in [T7]ll8j. 

In comparison with earlier studies, a more complete 
treatment of CM effects has been included in the model 
6, 29]. In the CM frame the transition operator becomes a 
two-body operator even in the case of a one-body nuclear 
current. As a consequence, the one-body current can give 
a contribution to the cross section of two-particle emis- 
sion independently of correlations. These effects were not 
properly taken into account in the previous calculations 
for proton-neutron knockout [9"[[T6]. Accounting for CM 
effects is not trivial since the lack of orthogonality between 
bound and scattering states (which are obtained from an 
energy-dependent optical potential) may give rise to spu- 
rious contributions to the calculated cross section. This 
issue has recently been overcome in [BJ enforcing orthog- 
onality between single particle initial and final states by 
means of the Gram-Schmidt procedure. The results show 
that the CM effects depend on kinematics. For the partic- 
ular case of super-parallel kinematics, which were used in 
the measurement of the 16 0(e,e'pn) 14 N reaction [9], they 
enhance the contribution to the cross section which arises 
from the one-body currents. This effect is dramatic at low 
missing momentum and fully accounts for the previously 
observed 9 discrepancy with experiment. The comparison 
between the 16 0(e,e'pn) 14 N data and the new calculations 
is reported in Sec. 14.11 



3 Experimental set-up 

3.1 The 16 0(e,e'pn) 14 N reaction 

A first measurement of 16 0(e,e'pn) 14 N reaction [5] was 
made at the electron scattering facility (3-spectrometer 
facility [30]) at MAMI, Mainz [3T[[32]. Data were taken 
with an incoming electron beam of energy 855 MeV at 
currents of 10-20 /iA. The beam was incident upon a wa- 
terfall target 33J of thickness 74 mg cm -2 . The data were 
collected at energy and momentum transfers of 215 MeV 
and 316 MeV/c where the ejected proton was detected in 
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the forward direction, parallel to q, with the ejected neu- 
tron detected in the backward direction, anti-parallel to q, 
in so called "super-parallel" kinematics. The ejected pro- 
ton and scattered electron were detected with Spectrom- 
eters A and B [3D] of the 3-spectrometer set-up while the 
ejected neutron was detected using the Glasgow-Tubingen 
time-of- flight detector system [34]. Further details about 
the experimental set-up and analysis of the data can be 
found in ref . [9] . The experimental resolution of the set-up 
was sufficient to distinguish groups of states in the resid- 
ual nucleus but not good enough to separate individual 
states. 



3.2 The 16 0(7,pn) 14 N reaction 

The 16 0(7,pn) 14 N reaction was measured at the Glasgow 
photon tagging facility [351136] at MAMI, Mainz [3TU32] . 
An electron beam of energy 855 MeV used at a current of 
50 nA was incident upon a 4 /im Nickel radiator to pro- 
duce tagged Bremsstrahlung photons in the energy range 
100 to 800 MeV. The Glasgow-tagger has an energy res- 
olution of 2 MeV. The tagged photons, collimated to a 
diameter of 18 mm, were incident upon a target of 1 mm 
thickness. The target cell was filled with deuterated water 
and consisted of an Aluminium frame with polythene foil 
windows of 30 /im thickness which was orientated at an 
angle of 30° with respect to the photon beam. 

The ejected protons were detected in an array of four 
hyper-pure Germanium detectors (HPGe) of the Edin- 
burgh Ge6- Array [37] , each of which covered a solid angle 
of 59 msr and had a proton energy acceptance of 18 - 
250 MeV. Pairs of double sided silicon strip detectors [21] 
positioned in front of the HPGe detectors were were used 
to determine the trajectory of the ejected protons and re- 
construct the reaction vertex. The ejected neutrons were 
detected at forward angles using the Glasgow-Tubingen 
time-of- flight detectors [34]. Five neutron detector stands 
were used which covered an in-plane polar angular range 
of 6 — 53° and a total solid angle of 146 msr. A pulse- 
height threshold of 5 MeV ee was used in the neutron de- 
tectors which resulted in a neutron kinetic energy thresh- 
old of «10 MeV. Full details of the experimental set-up 
and analysis of the data can be found in [21] . The exper- 
imental resolution of the set-up was not able to resolve 
individual excited states in the residual 14 N nucleus, only 
groups of states. 



4 Results 

4.1 The 16 0(e,e'pn) 14 N reaction 

Figure Q] shows the experimental and theoretical cross 
sections of the 16 0(e,e'pn) 14 N reaction as a function of 
the absolute magnitude of the missing momentum p m = 
q — p' p — p' n , where p' p and p' n are the momenta of the 
ejected nuclcons. The experimental cross section has been 
determined for a group of states in the residual 14 N for 
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Fig. 1. The 16 0(e,e'pn) 14 N cross section shown as a func- 
tion of the missing momentum for events in the range 2 < 
< 9 MeV for energy and momentum transfers of 215 MeV 
and 316 MeV/c. The curves show the results from theoretical 
calculations of the cross section which includes transitions to 
the first three excited states in 14 N, 2.31 MeV (0 + ), 3.95 MeV 
(1 + ) and 7.03 (2 + ). The dashed line is calculated only with the 
one-body currents; the dotted line also includes the 7r-seagull 
term; the dashed dotted includes the one-body, 7r-seagull term 
and pion-in-flight terms and the solid line is for the complete 
cross-section including contributions from IC. 

an excitation energy range of 2 to 9 MeV. The theoreti- 
cal curves are the result of DW calculations and are the 
average cross section of calculations for the kinematic set- 
tings as given in [9]. The calculations represent the sum 
of contributions for transitions to three excited states in 
14 N: the 2.31 MeV (0+), 3.95 MeV (1+) and 7.03 MeV 
(2+) states. 

The theoretical curves of fig. Q] also show the contri- 
butions of different terms of the nuclear current to the 
cross section. Cumulative contributions of the one-body, 
7r-seagull, pion-in-flight and isobar currents are all shown. 
At low missing momentum the largest contribution to 
the theoretical cross section is from one-body hadronic 
currents. Above p m = 150 MeV/c the 7r-seagull and ICs 
become increasingly more important with increasing p m . 
The pion-in-flight contribution is relatively small over the 
whole missing momentum range shown. 

The shape of the experimental and theoretical cross 
sections in fig. [1] show reasonable agreement in that they 
both decrease roughly exponentially with increasing p m 
and both show a flattening in the cross section at p m w 
175MeV/c. The magnitude of the two cross sections is 
in much better agreement compared to a previous com- 
parison in [9] where the theoretical calculations under- 
predicted the experimental data at lowp m . This improve- 
ment is due to the enhancement, at low p m , of the contri- 
bution from the one-body currents produced by the CM 
effects included in the present model [B]. 

Figure [2] shows a comparison of calculations of the full 
cross sections, including the one-body and two-body cur- 
rents, for transitions to the three different excited states 
included in the curves of fig. [1] The main strength in the 
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Fig. 2. Theoretical le O(e,e'pn) 14 N cross sections for energy 
and momentum transfers of 215 MeV and 316 MeV/c. The 
2.31 MeV (0 + ), 3.95 MeV (1 + ), 7.03 (2+) and the three states 
combined, represented by the dashed, dotted, dashed-dotted 
and solid lines respectively. The plots are for the full cross 
section including the one-body, 7r-seagull, pion-in-flight and IC 
terms. 

cross section is predicted to come from transitions to the 
3.95 MeV (1+) state up to p m w 290MeV/c where transi- 
tions to the 7.03 (2 + ) state become dominant. The calcu- 
lated contribution from transitions to the 2.31 MeV (0 + ) 
state is at least an order of magnitude weaker, over the 
full p m range shown, than those involving transitions to 
either of the other two states. 

The calculations in figs. 1 and 2 are performed in the 
DW approach for FSI. NN-FSI effects depend on kinemat- 
ics and on the reaction type and are generally small in 
proton-neutron emission [T7lll8j. For the 16 0(e,e'pn) 14 N 
reaction in the super-parallel kinematics NN-FSI are small 
but not negligible [6]. The effect of the mutual interac- 
tion between the two outgoing nucleons is shown in fig. [3l 
where the cross sections obtained in the DW and DW-NN 
approaches are compared for transitions to the 3.95 MeV 
(1 + ) state in 14 N. This one state dominates the reac- 
tion over nearly all of the measured p m range. The ef- 
fects of NN-FSI on the calculated cross section are rel- 
atively small. There is a slight decrease in cross section 
for p m < 50MeV/c and a slight increase for 150 < p m < 
225MeV/c and above p m = 300MeV/c. In general the 
calculations predict that NN-FSI have little importance 
for the kinematics shown here. This fact justifies the per- 
turbative treatment of NN-FSI. 



4.2 The 16 0(7,pn) 14 N reaction 

Figure S] shows the cross section for the 16 0(7,pn) 14 N as a 
function of the absolute magnitude of the missing momen- 
tum, p m , of the reaction. The data are shown for an inci- 
dent photon energy range of 150 < E~ < 250 MeV, proton 
in-plane azimuthal acceptance of 142 < 9 P < 158° and 
neutron in- plane azimuthal acceptance of 8 < 9 n < 32°. 
The experimental cross section has been determined for 
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Fig. 3. Theoretical 16 0(e,e'pn) 14 N cross sections for transi- 
tions to the 3.95 MeV (1 + ) excited state of 14 N for energy and 
momentum transfers of 215 MeV and 316 MeV/c. The solid 
curve uses the DW approach, the dashed line the DW-NN ap- 
proach for FSI. 
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Fig. 4. The 16 0(7,pn) 14 N cross section as a function of the 
missing momentum for events in the range 2 < < 10 MeV. 
The incident photon energy range was 150 < -E 7 < 250 MeV. 
The curves show the theoretical cross section for transitions to 
the 3.95 MeV (1 + ) state. The dashed line is calculated with 
only one-body currents included; the dotted line also includes 
the 7r-seagull term; the dashed dotted includes the one-body, 
7r-seagull term and pion-in-flight terms and the solid line is for 
the complete cross-section including contributions from IC. 



a group of states in the recoiling 1 N nucleus for an exci- 
tation energy range of 2 to 10 MeV. Figure 0] also shows 
the results of DW theoretical calculations for the reaction. 
The curves are for transitions to the 3.95 MeV (1 + ) state 
which is believed to dominate the cross section as in [14] 
and have been averaged over the kinematic settings which 
cover the acceptance of the experimental data. 

The theoretical curves of fig. 2] show the contributions 
of different terms of the nuclear current to the cross sec- 
tion. Cumulative contributions of the one-body, 7r-seagull, 
pion-in-flight and isobar currents are all shown, see the 
caption of fig. [4] for details. At low p m the largest contri- 



D.G. Middleton et al.: Knockout of proton-neutron pairs from O with electromagnetic probes 



5 




1(l -3 i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

lu 50 100 150 200 250 300 

p m [MeV/c] 

Fig. 5. Calculations of the le O(7pn) 14 N cross section for tran- 
sitions to the 3.95 MeV (1 + ) excited state of 14 N. The solid 
curve uses the DW approach, the dashed line the DW-NN ap- 
proach for FSI. 

bution to the theoretical cross section is from one-body 
hadronic currents. The inclusion of the 7r-seagull term 
causes a decrease in calculated cross section until roughly 
p m — 200 MeV/c where it has very little effect. The fur- 
ther inclusion of the pion-in-fiight contributions increases 
the cross section to roughly the same strength of the one- 
body hadronic current cross section for p m < 100 MeV/c 
after which point it increases the calculated cross section 
relative to the one-body hadronic currents alone. The in- 
clusion of ICs increases the calculated cross section for the 
whole p m range shown. 

Both the theoretical and experimental cross sections 
shown in fig. 4 show a similar trend of falling roughly ex- 
ponentially with increasing p m . The theory over-predicts 
the experimental data at low p m . Fig. 2] suggests that the 
discrepancy may decrease with increasing p m but more 
accurate measurements are necessary to confirm this. 

The effect of the mutual interaction between the two 
outgoing nucleons for the (7,pn) reaction is shown in fig. El 
Theoretical cross sections were obtained using the DW 
and DW-NN approaches for transitions to the 3.95 MeV 
(1+) state in 14 N. At low p m the effects of NN-FSI on 
the calculated cross section are very small. From about 
p m = 100 MeV/c the importance of NN-FSI increases un- 
til roughly p m — 200 MeV/c after which their importance 
again diminishes. This is in contrast to what was seen for 
the (e,e'pn) reaction where NN-FSI had very little effect 
on the calculated cross section. The inclusion of NN-FSI 
increases the theoretical cross section at high p m which, 
however, remains well within the statistical error bars as- 
sociated with the data points in this region. 



5 Conclusions 

An improved treatment of centre-of-mass effects in the 
electromagnetic current operator of the Pavia model has 
resulted in closer agreement with experimental data for 



the 16 0(e,e'pn) 14 N reaction with both the shape and mag- 
nitude of the experimental cross section being well de- 
scribed. 

A further comparison of the improved model with data 
from a recent measurement of the 16 0(7,pn) 14 N reaction 
showed a similar shape with p m but over-predicted the 
strength of the measured cross section at lowp m . However, 
more accurate data would be required in order to draw a 
more detailed interpretation from this comparison. 
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